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ABSTRACT 

Aims. We introduce our imaging survey of possible young massive globular clusters in M31 performed with the Wide Field and 
Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). We obtained shallow (to B~ 25) photometry of individual stars 
in 20 candidate clusters. We present here details of the data reduction pipeline that is being applied to all the survey data and describe 
its application to the brightest among our targets, van den Bergh (VdBO), taken as a test case. 

Methods. Point spread function fitting photometry of individual stars was obtained for all the WFPC2 images of VdBO and the 
completeness of the final samples was estimated using an extensive set of artificial stars experiments. The reddening, the age and the 
metallicity of the cluster were estimated by comparing the observed color magnitude diagram (CMD) with theoretical isochrones. 
Structural parameters were obtained from model-fitting to the intensity profiles measured within circular apertures on the WFPC2 
images. 

Results. Under the most conservative assumptions, the stellar mass of VdBO is M > 2.4 x 10 4 Mq, but our best estimates lie in the 
range ^ 4 - 9 x 10 4 Mq. The CMD of VdBO is best reproduced by models having solar metallicity and age = 25 Myr. Ages less than 
^12 Myr and greater than 60 Myr are clearly ruled out by the available data. The cluster has a remarkable number of red super 
giants (^ 18) and a CMD very similar to Large Magellanic Cloud clusters usually classified as young globulars such as NGC 1850, 
for example. 

Conclusions. VdBO is significantly brighter (^ 1 mag) than Galactic open clusters of similar age. Its present-day mass and half-light 
radius (r h = 7.4 pc) are more typical of faint globular clusters than of open clusters. However, given its position within the disk of 
M31, it is expected to be destroyed by dynamical effects, in particular by encounters with giant molecular clouds, within the next ~ 4 
Gyr. 

Key words. Galaxies: star clusters - Galaxies: individual: M31 - (Stars:) supergiants - Stars: evolution 
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1. Introduction 

Much of the star formation in the Milky Way is thought to have 
occurred within star clusters (Lada et al. |199U Carpenter et al. 
2000). Therefore, understanding the formation and evolution of 
star clusters is an important piece of the galaxy formation puzzle. 
Our understanding of the star cluster systems of spiral galaxies 
has largely come from studies of the Milky Way. Star clusters 
in our Galaxy have traditionally been separated into two vari- 
eties, open and globular clusters (OCs and GCs hereafter). OCs 
are conventionally regarded as young (< 10 10 Gyr), low-mass 
(< 10 4 Mq) and metal-rich systems that reside in the Galactic 
disk. In contrast, GCs are characterized as old, massive systems. 
In the Milky Way, GCs can be broadly separated into two com- 
ponents: a metal-rich disk/bulge subpopulation, and a spatially 
extended, metal -poor halo subsystem (Kinman l 19591 Zinn[l985; 
see also Brodie & Strader 20061 Harris[2001 for general reviews 
of GCs). 

However, the distinction between OCs and GCs has become 
increasingly blurred. For example, some OCs are sufficiently lu- 
minous and old to be confused with GCs (e.g., Phelps & Schick 
2003). Similarly, some GCs are very low-luminosity systems 
(e.g., Koposov et al. 120071 1 and at least one has an age that is 
consistent with the OC age distribution (Palomar 1; Sarajedini 
et al. 120071 1. Moreover, a third category of star cluster, "young 
massive clusters" (YMCs) are observed to exist in both merg- 
ing (e.g., Whitmore & Schweizer 1 19951 1 and quiescent galax- 
ies (Larsen & Richtler 1999), Indeed, YMCs have been known 
to exist in the Large Magellanic Cloud for over half a century 
(Hodge |19611 l. These objects are significantly more luminous 

than OCs (My -8 up to My 15), making them promising 

candidate young GCs. Once thought to be absent in the Milky 
Way, recent observations suggest that their census may be quite 
incomplete, as some prominent cases have been found recently 
in the Galaxy as well (Clark et al. [2005; Figer |2008l >. 

Thus, a picture has emerged that, rather than representing 
distinct entities, OCs, YMCs and GCs may represent regions 
within a continuum of cluster properties dependent upon local 
galaxy conditions (Larsen 120031 1. The lifetime of a star cluster 
is dependent upon its mass and environment. Most low-mass 
star clusters in disks are rapidly disrupted via interactions with 
giant molecular clouds (Lamers & Gieles 2006; Gieles et al. 
2007). These disrupted star clusters are thought to be the ori- 
gin of much of the present field star populations (Lada & Lada 
120031) . Surviving disk clusters may then be regarded as OCs or 
YMCs, depending upon their mass. Star clusters in the halo may 
survive longer since they are subjected to the more gradual dy- 
namical processes of two-body relaxation and evaporation. The 
clusters which survive for an Hubble time - more likely to occur 
away from the disk - are termed GCs (see also Krienke & Hodge 
120071 1. To date, no known thin disk GCs have been identified in 
the Milky Way. 

After the Milky Way, M31 is the prime target for expand- 
ing our knowledge of cluster systems in spirals. However, our 
present state of knowledge about the M3 1 cluster system is far 
from complete. Similar to the Milky Way, M31 appears to have 
at least two GC subpopulations; a metal-rich, spatially concen- 
trated subpopulation of GCs and a more metal-poor, spatially 
extended GC subpopulation (Huchra et al. 1991 ; Barmby et al. 
2000). Also, again similar to the Milky Way GCs, the metal-rich 
GCs in M3 1 rotate and show "bulge-like" kinematics (Perrett et 
al. |2002| l. However, unlike the case in the Milky Way, the metal- 
poor GCs also show significant rotation (Huchra et al. 1 199 1 1 
Pen-ett et al. |2002l Lee et al. |2008l . Using the Perrett et al. ( [20021 



data, Morrison et al. (2004) identified what appeared to be a thin 
disk population of GCs, constituting some 27% of the Perrett et 
al. (|2002) sample. Subsequently, it has been shown that at least 
a subset of these objects are in fact young (< 1 Gyr), metal-rich 
star clusters rather than old "classical" GCs (Beasley et al. 2004; 
Burstein et al. 120041 Fusi Pecci et al. l2005l Puzia et al. I2005I 

Fusi Pecci et al. ( 120051 hereafter F05) presented a compre- 
hensive study of bright young disk clusters in M31, selected 
from the Revised Bologna Catalogue^ (RBC, Galleti et al. |2004l 
by color [(B - V) < 0.45] or by the strength of the HB line in 
their spectra (HB > 3.5A). While these clusters have been noted 
since Vetesnik ( 1962) and have been studied by various authors, 
a systematic study was lacking. F05 found that these clusters, 
that they termed - to add to the growing menagerie of star clus- 
ter species - "Blue Luminous Compact Clusters" (BLCCs), are 
fairly numerous in M31 (15% of the whole GC sample), they 
have positions and kinematics typical of thin disk objects, and 
their colors and spectra strongly suggest that they have ages (sig- 
nificantly) less than 2 Gyr. 

Since they are quite bright (-6.5 S My $ -10.0) and - at 
least in some cases - morphologically similar to old GCs (see 
Williams & Hodge |200T| hereafter WH01), BLCCs could be re- 
garded as YMCs, that is to say, candidate young globular clus- 
ters. In particular, F05 concluded that if most of the BLCCs have 
an age ^ 50 - 100 Myr they are likely brighter than Galactic 
Open Clusters (OC) of similar ages, thus they should belong to 
a class of objects that is not present, in large numbers, in our 
own Galaxy. Unfortunately, the accuracy in the age estimates 
obtained from the integrated properties of the clusters is not suf- 
ficient to determine their actual nature on an individual basis, 
i.e., to compare their total luminosity with the luminosity dis- 
tribution of OCs of similar age (see Bellazzini et al. 2008] and 
references therein). 

In addition to the question of the masses and ages of these 
BLCCs, it has become clear that the BLCC photometric and 
spectroscopic samples in M3 1 may suffer from significant con- 
tamination. Cohen, Matthews & Cameron (2006, hereafter C06) 
presented NIRC2@KeckII Laser Guide Star Adaptive Optics 
(LGSAO) images of six candidate BLCCs. Their K' very-high 
spatial resolution images revealed that in the fields of four of the 
candidates there was no apparent cluster. This lead C06 to the 
conclusion that some/many of the claimed BLCC may in fact be 
just asterisms, i.e. chance groupings of stars in the dense disk 
of M3 1 . While the use of the near infrared K' band (required by 
the LGSAO technique) may be largely insensitive to very young 
clusters that are dominated by relatively few hot stars, which 
emit most of the light in the blue region of the spectrum, the in- 
ference is that the true number of massive young clusters of M3 1 
may have been severely overestimated. 

Therefore, in order to ascertain the real nature of these 
BLCCs we have performed an HST survey to image 20 BLCCs 
in the disk of M31 (program GO-10818, PL: J. Cohen). The key 
aims of the survey are: 

1 . to check if the imaged targets are real clusters or asterisms, 
and to determine the fraction of contamination of BLCCs by 
asterisms; 

2. to obtain an estimate of the age of each cluster in order to 
verify whether it is brighter than Galactic OCs of similar age. 
Ultimately the survey aims to provide firm conclusions on 
the existence of BLCCs (YMCs) in M3 1 as a distinct class 



1 www.bo.astro.it/M31 
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Fig. 1. F450W mosaic of the whole field sampled by our WFPC2 
observations. The cluster VdBO is at the center of the PC camera. 



of object with respect to OCs (see Krienke & Hodge 120071 
120081 and references therein). 

In the present contribution we describe the data reduction 
and analysis strategies that we will apply to our cluster sample to 
estimate their ages and metallicities. The overall procedure is de- 
scribed using the brightest among the observed clusters, VdBO, 
as a specific case. We conclude this section with a brief presen- 
tation of the cluster VdBO, below. 

The present paper is organized as follows. The observa- 
tions and the data reduction procedure are described in detail 
in Sect. 2; the principal assumptions that will be adopted in the 
whole survey are also reported in this section. Sect. 3 is devoted 
to the analysis of the surface brightness profile and of the Color 
Magnitude Diagram of VdBO, including total luminosity, age 
and metallicity estimates. In Sect. 4 our main results are briefly 
summarized and discussed. 

1.1. The cluster van den Bergh (VdBO) 

VdBO was indicated as an open cluster by Hubble d!936l l in 
the image on the frontispiece of his book The Realm of the 
Nebulae^ van den Bergh (|1969l l presents VdBO as the bright- 
est open cluster of M31, reporting an integrated spectral type 
AO. He also notes that the cluster contains the Cepheid variable 
V40 (Hubble [19291 . A check of Hubble's ( [19291 finding charts 
revealed that two sources are labeled # 40 in his plate VII: one of 
them seems indeed associated with the cluster, while the other is 
~ 8' away from VdBO , near the association OB78 = NGC 206 
(van den Bergh [T9641 see also Hodge [T979I I. The cluster was 
re-discovered by Hodge 1979, who classified it as an open clus- 
ter (C107, see also HodgeQMQ. Finally, Battistini et al. ( fl987T > 
listed the cluster as their class D candidate globular cluster num- 
ber 195 (B195D in the RBC). The failure to identify B195D with 

2 S. van den Bergh kindly drove our attention to this curious occur- 
rence. 



VdBO was due to the fact that the coordinates provided by van 
den Bergh ( 1 19691 ) were in error by 17". For this reason VdBO 
and B195D survived as independent entries in M31 GC cata- 
logues until the present day. In our survey we imaged both the 
clusters and the WFPC2 images revealed unequivocally that the 
two targets are in fact the same cluster. In particular the images 
intended to observe B195D have the cluster in the center of the 
PC camera while in the VdBO images the cluster lie in the cor- 
ner of the PC opposite to the WF cameras, such that part of the 
cluster is out of the image. In the following (and in the future) 
we will refer to the cluster as VdBO. The dataset analysed here 
is the one with the cluster centered on the PC images, hence the 
actual label in the header of the fits files is B195D. 

VdBO is located at a projected distance of R p = 10.8 kpc 
from the center of M31 to the South-West, just ~ 4' from the 
major axis of the galaxy (see Tab. [1), near the edge of one of 
the most prominent substructures of the M3 1 disk, the so called 
10 kpc ring (see Hodge 1992 and Barmby et al. 2006 and refer- 
ences therein) and within a the large OB association OB80 (van 
den Bergh 1 19641 A80 in Hodge (fTMll l atlas). Its radial velocity 
(V r = -567 km/s, Perrett et al. 2002) is in full agreement with 
the rotation curve of the HI disk of M31 (Carignan et al. 2006), 
thus confirming the physical association with the thin disk of 
the parent galaxy (F05). The strong value of the Hp index sup- 
ports the idea that the cluster is younger than 1 Gyr (Hp = 4.3 
A, Perrett et al. [20020). The existing estimates of both V r and 
Hp are nicely confirmed by recent high signal-to-noise spectra 
acquired at the Italian Telescopio Nazionale Galileo (S. Galleti, 
private communication). 

With the assumed reddening and distance, the integrated V 
magnitude reported in the RBC (see Tab. [T]( gives an absolute 
magnitude My = -10.03, much brighter than any Galactic open 
cluster older than 10 Myr (see Bellazzini et al. 120081 and below); 
it appears quite extended and irregular in shape even in ground 
based images. In these ways VdBO stands out among the mem- 
bers of our candidate BLCC sample that are, in general, fainter 
and more compact than it. 

2. Observations and Data Reduction 

Our survey was originally planned for the Advanced Camera 
for Surveys (ACS) but it was performed with the Wide Field 
and Planetary Camera 2 (WFPC2) during cycle 16 because 
of the failure of ACS. For each target of our survey we ac- 
quired two F450W and two F814W images, all with 400 s ex- 
posure time and gain = 7e~/DN. The pointings were chosen to 
place the main target at the center of the PC (800 x 800 px 2 , 
with pixel scale 0.045 arcsec/px), while the three WF cameras 
(800 x 800 px 2 , with 0.099 arcsec/px) are supposed to sample 
the surrounding fields. The images of VdBO discussed here were 
acquired on July 2, 2007. The image of the whole WFPC2 mo- 
saic image is shown in Fig. 1 . It is clear that there are substruc- 
tures and density gradients on the scale of the whole mosaic im- 
age, mainly due to the inclusion of the edges of the large stellar 
association embedding the cluster (A80, Hodge 1 1 98 It - As the 
overall stellar density on the WF2 field is larger than in WF3 
and WF4, we make the conservative choice to adopt the WF2 
as our preferred sample of the background population that is 

3 Note that Perret's et al. measures refers to B195D, i.e. the "alter 
ego" of VdBO whose available coordinates were the most appropriate 
for the cluster. In this context, it is interesting to note that, adopting 
a calibration based on old GCs, Perrett et al. found [Fe/H]=-1.64 for 
VdBO, from integrated spectral indices (see F05). 
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Table 1. Positional and Photometric parameters for VdBO from the RBC" 



NAME 


alt NAME 


RA/2000 


Decj2ooo 


X 


Y 


U 


B 


V 


R 


J 


H 


K 


VdBO 


B195D" 


00:40:29.3 


+40:36:14.7 


-47.2' 


-4.3' 


14.97 


15.31 


15.06 


14.92 


13.77 


13.14 


12.99 



a X and Y are projected coordinates in the direction along (increasing Eastward) and perpendicular to the major axis of M31 (increasing 

Northward) respectively, in arcmin, see Galled et al. 2004, and references therein. 
b see Sect. 2.5. 



expected to contaminate the Color Magnitude Diagram of the 
cluster, while we will consider the average density over all the 
WF fields when we will compute stellar density profiles based 
on star counts (Sect. 3). In the present context, when we speak 
of "background population" we refer to all the stars belonging 
to the field of M3 1 but unrelated to the cluster we are studying. 
Zoomed views of the PC field in both F450W and F814W pass- 
bands are shown in Fig. [2] 

As the observational material and the degree of crowding are 
essentially the same for all the surveyed fields, we tuned our 
data-reduction strategy to be exactly the same in all cases, to 
maintain the highest degree of homogeneity in the final products 
of the survey. Data reduction has been performed on the pre- 
reduced images provided by STScI, using HSTPHOlQ(Dolph in 
2000a), a Point Spread Function -fitting package specifically de- 
voted to the photometry of WFPC2 data. The package identifies 
the sources above a fixed flux threshold on a stacked image and 
performs photometry on individual frames, and automatically 
applies the correction for the Charge Transfer Efficiency (CTE, 
Dolphin 2000b). It then transforms instrumental magnitude to 
the VEG AMAG system (see Holtzman et al. ( 119951 ) and Dolphin 
(2000b)), deals with cosmic-ray hits, and takes also into account 
all the information about image defects that is attached to the 
observational material. We fixed the threshold for the search of 
sources on the images at 3 cr above the background. HSTPHOT 
provides as output the magnitudes and positions of the detected 
sources, as well as a number of quality parameters for a suitable 
sample selection, in view of the actual scientific objective one 
has in mind. Here we selected all the sources having valid mag- 
nitude measurements in both passbands, global quality flag = 1 
(i.e., best measured stars), crowding parameter < 0.3, x 1 < 2.0 
and -0.5 < sharp < 0.5, in both passbands, (see Dolphin 2000a 
for details on the parameters). This selection cleans the sam- 
ple from the vast majority of spurious and/or badly measured 
sources without significant loss of information, and it has been 
found to be appropriate for the whole survey. 

In Fig.|3]the Color Magnitude Diagrams (CMD) of the fields 
imaged by the four chips of WFPC2 are shown. The threshold 
for the saturation of bright stars and the boundaries at which the 
completeness of the sample reaches 90%, 70% and 50% are also 
shown, as derived from the artificial stars experiments described 
below. As the CMD is quite typical of our survey, it is worthy of 
some general comments while a detailed analysis is deferred to 
Sect. 3 below. First, our photometry is relatively shallow, due the 
short exposure times of our images; the 50% completeness level 
is reached at F450W - 25.50. For the same reason our images, 
and particularly the F450W ones in which the background light 
is very low, are badly affected by CTE (see Fig.|2]i. Therefore the 
accuracy of the absolute and relative photometry is not particu- 



See http : //purcell . as . arizona . edu/hstphot/ 
5 Except for the very crowded region at the center of the cluster. For 
10 px < r < 50 px, the 50% completeness level is reached at F450W <: 
23.5. 



Table 2. Uncertainties in the relative photometry from artificial 
stars experiments, for 10 px < r < 160 px, PC field. 



F450W 


cr" 


F814W 


cr" 


18.00 


0.009 


18.00 


0.010 


18.50 


0.010 


18.50 


0.011 


19.00 


0.010 


19.00 


0.012 


19.50 


0.011 


19.50 


0.013 


20.00 


0.013 


20.00 


0.016 


20.50 


0.016 


20.50 


0.020 


21.00 


0.018 


21.00 


0.026 


21.50 


0.023 


21.50 


0.036 


22.00 


0.029 


22.00 


0.050 


22.50 


0.039 


22.50 


0.068 


23.00 


0.054 


23.00 


0.087 


23.50 


0.076 


23.50 


0.138 


24.00 


0.107 


24.00 


0.218 


24.50 


0.153 


24.50 


0.336 


25.00 


0.241 


25.00 


0.377 


25.50 


0.309 


25.50 


0.400 



cr are ± 1 standard deviations after the clipping of outliers at more 
than 3<r from the mean. 



larly good (see, for example, Fig.|4]and Tab.|2l below). In spite of 
that, the very wide wavelength baseline provided by the F450W 
and F814W filters produces relatively well defined sequences in 
the CMD (compare, for example, with the CMD of similar fields 
obtained by WH01 with the same camera and longer exposure 
times but using F439W and F555W filters). 

All the fields targeted by our survey cross the outer regions 
of the star-forming thin disk of M3 1 (see F05), and as a conse- 
quence, in most cases, the most prominent feature of the CMD 
is the nearly vertical plume of young Main Sequence stars that 
is seen in Fig.[3]around F450W - F8UW - 0.2. The wide blob 
of stars at F450W > 24.0 and F450W - F814W > 1.5 is con- 
sistent with being due to the brightest Red Giants near the tip of 
the Red Giant Branch (RGB) of the old-intermediate population 
that seems to be pervasive in the M31 disk (see Bellazzini et al. 
2003, and references therein). Red and blue supergiants as well 
as other less-massive evolved stars are likely present at bright 
magnitudes over the whole color range covered by our CMD 
(see Massey |2006l) . 

2. 1 . Artificial stars experiments 

The completeness of the samples and the accuracy in the relative 
photometry are best estimated with extensive sets of artificial 
stars experiments (see Bellazzini et al. 2002a 2002b and Tosi et 
al. l2001l for detailed discussions and references). 

HSTPHOT allows easy, fast and fully automated runs of arti- 
ficial stars experiments. Fake stars in a user-selected color range, 
extracted at random from a Luminosity Function (LF) similar to 
the observed one, are added to the original frames one at a time 
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Fig. 3. CMD of the fields sampled by the four chips of the WFPC2. The number of stars plotted is reported in the upper left corner 
of each panel. The upper line marks the threshold above which stars saturate the intensity scale of the images. The lower lines are 
CMD loci at the same level of completeness, 90%, 70% and 50% from top to bottom, respectively (see labels in the WF2 panel). 



to avoid self-crowding (Dolphin, private communication) and 
the photometric reduction is repeated. With the final catalogue 
of input and output magnitudes of artificial stars the distribution 
of photometric errors and the completeness of the samples can 
be studied as a function of color and as a function of the distance 
from the center of the cluster under consideration (i.e. as a func- 
tion of crowding). We simulated a total of 728398 artificial stars, 
roughly equally distributed on the four WFPC2 chips. 

Fig. |4] shows the distributions of the differences between 
the output and input magnitudes of artificial stars as a func- 
tion of F450W (left panels) and F814W (right panels) magni- 
tudes, providing a direct estimate of the typical uncertainties of 
our relative photometry. The small excess of stars at negative 
m ou r - ffii n , increasing in number and amplitude of the differ- 
ence for fainter magnitudes, is due to artificial sources that are 
erroneously recovered with a brighter magnitude because they 
are blended with real sources present on the image (see Tosi 



et al. 120011) . Even in the most crowded region of the PC that 
includes the cluster (top panels of Fig. HJ the effects of blend- 
ing are not particularly severe, at least for relatively bright stars. 
The probability of a star with F450W < 23.5 to have its mag- 
nitude decreased by more than 0. 1(0.2) mag by the combination 
of blending and photometric error is 2.8%(1.4%) if its color lies 
in the range -0.6 < F450W - F814W < 1.5 and 3.5%(1.6%) 
for 2.0 < F450W - F8UW < 4.0. Typical photometric uncer- 
tainties as a function of magnitude are reported in Table [2] for 
the innermost region of the PC field, covering most of the cluster 
that is the main subject of the present study. 

Finally the completeness factors (C/) as a function of mag- 
nitude for different regions of the PC and for the WF fields 
are shown in Fig. [5] for stars in the wide color range -0.6 < 
F45QW - F814W < 1.5. Outside of the innermost region of the 
PC including the cluster, the Cf functions are nearly indistin- 
guishable. For r > 50 px the completeness is larger than 80% 
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Fig. 4. Distributions of the differences between the output and input magnitudes of artificial stars as a function of F450W (left 
panels) and F814W (right panels) magnitudes, for the PC and WF fields. The top panel displays the distributions for the most 
crowded region of the PC camera, i.e. the one containing the cluster, r^j is the distance from the cluster center in PC pixel units, 
assuming (x,y)=(405,398) as the coordinate of the center in the reference frame of the photometric catalogue. To make the diagrams 
more easily readable we plot just a fraction of the whole set of artificial stars, i.e. 50000 stars per field, approximately, while more 
than 150000 per field are typically recovered. 



for f450W < 24.0 and in any case Cf - 1 (i.e. completeness 
100%) for F450W < 22.0. 

2.2. Theoretical stellar models 

Most of our inferences about the physical parameters of the stel- 
lar populations (clusters or field) considered in our survey will be 
obtained from the comparison between the observed CMDs and 
theoretical stellar models, in the form of isochrones or synthetic 
CMDs. The need to have models in the natural photometric 
system in which the observations were obtained (HST/WFPC2 
VEGAMAG) and to have a set of isochrones reaching ages as 
young as 10 Myr led us to chose the set by Girardi et al. (2002, 
hereafter G02), as our reference grid of stellar models. In partic- 



ular we took their HST-color version of the solar- scaled models 
by Salasnich et al. (2000), with overshooting and a simplified 
TP-AGB evolution, as this set includes 10 Myr old isochrones 
up to super-solar metallicitiefl In some cases, when a particu- 
lar model is needed, we use the CMD web tooQ (Marigo et al. 
2008 ), that allows the on-line computation of models from user 
specified inputs, using the G02 set. 

In some cases, for comparison and/or for special applica- 
tions, we use the BAST0 database, collecting the theoretical 
models by Pietrinferni et al. ( 120041 ). and updates. In particu- 



6 http://pleiadi.oapd.inaf.it 

7 http : // stev . oapd . inaf . it~lgirardi/ cgi-bin/ cmd 

8 http://www.oa-teramo.inaf.it/BASTI/index.php 
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Fig. 2. F450W (upper panel) and F8 1 4W (lower panel) images of 
the whole PC camera, with VdBO at the center. The superposed 
circles have radius r=160, 205, 260, 288 and 330 pixels, from 
inside out, and mark the edges of the annuli whose CMDs are 
shown in Fig. [9] below. The light stripes associated with stars in 
the F450W image are due to the effect of CTE that is particularly 
strong in this shallow low-background image. 



lar BASTI provides a very practical Web Tool to produce syn- 
thetic CMDs of populations with ages, chemical composition, 
initial mass function, binary fraction (/j) etc. selected by the 
user (Cordier et al. 2008), that can be used to compare mod- 
els and observations in term of star counts in different color and 
magnitude ranges (see Fig. [6j for an example of application). 
Unfortunately, the models are not provided in the WFPC2 pho- 



KKI4M .. Fig- 5. Completeness factor (Cf) as a function of F450W mag- 
nitude for the color range enclosing MS stars. Upper panel: C/ 
for regions of the PC field at different distances from the cluster 
center. Lower panel: Cf for the three WF fields. Note that the 
three curves are indistinguishable within the uncertainties. 



tometric system - so theoretical magnitudes have to be trans- 
formed - and isochrones/synthetic CMDs for ages < 30 Myr are 
not provided; for these reasons we didn't adopt the BASTI set as 
the reference for our survey. In the considered range of ages G02 
and BASTI isochrones (with overshooting) provide very similar 
predictions of color and magnitudes, while evolving masses may 
differ by ~ 20% (see also Gallart, Zoccali & Aparicio 2005). 



2.3. Reddening and Distance 

To correct for the effects of interstellar extinction and reddening 
we will always adopt the relations A F ^ W = 4.0 \5E(B - V) and 
Af 8 i4 = 1.94 8ff(l? - V), as reported by Schlegel, Finkbeiner & 
Davis d 1 9981 >. As our clusters are embedded in the structured 
dusty disk of M31 it does not seem appropriate to assume a 
unique value of reddening for all of them; the typical redden- 
ing value attributed to Galactic dust toward M31 ranges from 
E(B - V) = 0.06 (Schlegel et al. [19981 to E(B - V) =* 0.11 
(see Galleti 120041 and references therein), but it is likely that 
our clusters are more reddened than this (Barmby et al. 2000; 
Fan et al. 2008). To get an estimate of the reddening affect- 
ing the clusters in our survey we compare theoretical models 
(isochrones and synthetic CMDs) to the observed MS in the 
range 22.0 S5 F45QW S5 24.0. In this range, corresponding to 
absolute magnitudes -3.0 5 Mf45ow ^ 0.0, the color of the MS 
is only weakly sensitive to metallicity and various sets of theoret- 
ical models provide very consistent predictions. An example of 
our analysis is presented in Fig. [6] where we compare the color 
distribution at the blue edge of the MS of the observed sample 
and of synthetic samples (from the BASTI webtool) of different 
metallicities, adopting different reddening values. The compar- 
isons confirm that the sensitivity to metallicity of the reddening 
estimate is very weak, as expected. In the case of VdBO we ob- 
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Fig. 6. The observed CMD of VdBO (black dots, only stars with 
r < 160 px) is compared with the synthetic CMD (grey x sym- 
bols) of 30 Myr old, f b =50% populations having Z = 0.008 
(left panel), Z = 0.019 (middle panel), a ndZ = 0.040, obtained 
from the BASTI webtool (Cordier et al. 2008), transformed to 
WFPC2-VEGAMAG with Dolphin (I2000bb equations, and cor- 
rected for photometric errors and completeness according to the 
results of our artificial stars experiments. The thin lines enclose 
the selection box in which the cumulative color distributions 
shown in the upper panels have been obtained, focusing on the 
blue edge of the Main Sequence. In these panels the observed 
color distribution (continuous line) is compared to the distribu- 
tions of the synthetic sample of the adopted metallicity for three 
different assumptions on the reddening value (dashed lines), re- 
ported in the upper label. The middle value corresponds to the 
distribution that best fits the observations and is also reported in 
the upper left corner of the CMDs. Note the very weak depen- 
dence of the reddening estimate on the metallicity of the adopted 
model. 



Fig. 7. Isochrones of different ages and metal content are plot- 
ted on the "visibility window" of our CMDs, enclosed on the 
bright side by the saturation limits and on the faint side by the 
Cf = 50% line (long-dashed lines). The continuous curves are 
isochrones from the G02 set; ages and metallicities are indicated 
in the figure. 



their Turn Off (TO) points and/or from the distribution of their 
Super Giant populations. 

In the metallicity range that is most likely to enclose the disk 
populations (we are considering |Zq $ Z 5 2Zq) we can detect 
the TO point of clusters roughly ranging from 10 to 500 Myr old. 
As the only BLCCs for which a direct CMD-based age estimate 
has been obtained are 60-160 Myr old (WH01), the age sensi- 
tivity of the survey seems rather appropriate; however clusters 
in the age range 0.5 - 2 Gyr may prove very difficult to age date 
with our data. For the oldest populations (age ^ 2 Gyr) we can 
hope to detect just the tip of the RGB, as shown by the age=12 
Gyr isochrones plotted as thick lines in Fig. [9] below. 



tain E(B — V) = 0.2 ± 0.03 with this method, and we will always 
adopt this value below. 

In the following and for the whole survey we adopt (m - 
M)q = 24.47 ± 0.07 as the distance modulus of all the consid- 
ered populations, from McConnachie et al. (2005), correspond- 
ing to an heliocentric distance D = 783 kpc. At this distance 1" 
corresponds to 3.8 pc, V to 228 pc. 

2.4. Accessible age range 

As the degree of crowding of all the surveyed fields is quite sim- 
ilar and the observational set-up is identical in all cases, the sat- 
uration limit and the Cf - 0.50 limit reported in the CMDs of 
Fig EJcan be considered representative of the typical CMD win- 
dow that is accessible with the survey data. In Fig.Qwe compare 
isochrones of different ages and metallicities with this window to 
have an idea of the age range in which we can obtain reasonable 
age estimates for the considered clusters from the luminosity of 



3. The CMD and structure of the cluster VdBO 

3. 1 . Distribution of resolved stars 

To identify the stellar population of the cluster as securely as 
possible, it is useful to have an idea of the surface density distri- 
bution of its resolved stars. In the present context we are inter- 
ested only in defining the characteristic size of the region domi- 
nated by cluster stars, in order to select samples of likely cluster 
members by radius (see Sect. 3.4 for a detailed analysis of the 
light profiles). 

Stars were selected on the CMD from the box shown in the 
diagram enclosed in the lower panel of Fig. [8] The box is ex- 
pected to pick up the best-measured MS and SG stars typical 
of the cluster population, while excluding populations that are 
clearly not associated with the cluster, such as the much older 
stars around the tip of the RGB. For r < 3" star counts are signif- 
icantly affected by radially varying incompleteness in the range 
of magnitudes considered. Beyond this limit the degree of com- 
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Fig. 8. Upper panel: Background-subtracted surface density pro- 
file of VdBO computed by counting stars on circular concentric 
annuli around the center of the cluster. The arrow marks the ra- 
dius where a sudden change of slope in the profile appears, at 
r ^ 160 px = 7.3". Lower panel: Background-subtracted pro- 
file from star-counts (filled circles with errorbars) converted to 
a major-axis profile, adopting the reported values of PA and e. 
Open symbols are the corresponding light profiles described in 
Sect. 3.4, squares for the F450W image and triangles for F8 14W, 
vertically shifted by an arbitrary normalization to match the star 
counts at r ma > 3". The dotted lines mark the average surface 
density in each of the WF cameras, the dashed line is the av- 
erage of the three, which was in the end adopted as the back- 
ground value to subtract to star-count profiles. Only stars within 
the L-shaped box plotted in the CMD in the upper right corner of 
the lower panel are selected for star counts, as probable cluster 
members. 

pleteness is fairly high and essentially constant with radius (see 
Fig- El above), hence the derived profile should be reliable. 

In the upper panel of Fig. [8] we show the surface density 
profile obtained by counting stars on circular annuli centered 
on the cluster center. The observed profile displays an obvious 
break at r 7.3", where it begins to decline with a gentler 
slope out to r ~ 14". The break in the profile may reflect an 
inner core + outer corona structure of VdBO, which is typical 
of Galactic Open Clusters (see Kubiak et al. 119921 Kharchenko 
et al. 120051 Mackey & Gilmore|2003 Elson et al. [19871 and ref- 
erences therein), or it may be — at least partially- due to the 
elongated distribution of the cluster stars unaccounted for by 
our adoption of circular annuli. To investigate this possibility we 
transformed the radial coordinate of each star (r) into a major- 
axis radius (r ma ) defined as 



= A PS + 



1 



(1) 



r 1(1 ~e) 
where 

X r = (X- X )cos(PA XJ ) + (Y - Y )sin(PA XJ ) (2) 



and (Xr,, Yq) are the coordinate of the center of the cluster, e = 
1 - b/a, is the ellipticity, where a and b are the semi-major and 
semi-minor axis, respectively, and PA X j is the position angle 
measured from the X axis toward the Y axis. Both e and PA X y 
are taken (or easily derived, in the case of PA X j) from the results 
of the analysis of the light distribution presented in Sect. 3.4, 
below. Eq.Q]has been adapted to our case from Eq. 4 by Martin 
et al. (120581 . 

The ellipticity-corrected major axis profile is plotted in the 
lower panel of Fig. [8] and it clearly shows that the change of 
slope in the original profile was an artifact due to the inadequacy 
of the assumption of circular symmetry. The result is supported 
by the good match between the star-counts profile and the light 
profiles (from Sect. 3.4) over the large radial range where they 
can be compared (r > 3"). 

It is interesting to note that the cluster profile appears to 
extend to remarkably large distances from the center, out to 
- 15" - 57 pc. As the process of profile analysis described in 
Sect. 3.4 includes also the fitting of King ( 1966 hereafter K66) 
models, it is interesting to note that the limiting radius of the 
K66 models that best fits the surface brightness profiles is also 
r, ^ 15", thus supporting the conclusion that the cluster is very 
extended. 

The elongated shape of the cluster will be taken into account 
in the detailed analysis of the profiles of Sect. 3.4. For present 
purposes it is sufficient to conclude that most of the cluster stars 
are enclosed within a (circular) radius of 7.3" (160 px) from the 
center. We take this as a reference radius for the following anal- 
ysis of the CMD, as it allows a very simple radial selection, re- 
membering that some cluster members are also present at larger 
radii. 

The upper left panel of Fig. |9]shows the CMD of stars within 
10 < r < 160 px, an annulus that, as stated earlier, should be 
dominated by cluster stars. The innermost r < 10 px region 
has been excluded because of severe incompleteness. A main se- 
quence with a TO around F450W ~ 21.5 is the most populated 
branch of the diagram, with a blue edge at F4-50W - F814W - 
0.0. Blue and red supergiants (BSGs, RSGs) are clearly identi- 
fied, spanning a large color range (0.0 < F4-5QW - F814W < 
3.6 mag). A 25 Myr isochrone of solar metallicity (from the 
G02 set) seems to provide a satisfactory fit to the MS and to the 
sizable luminosity range spanned by supergiants, suggesting an 
extended Blue Loop phase (see Williams & Hodge 2001). The 
color of the reddest supergiants is not fully reproduced (a long 
standing and not-so-critical problem of theoretical models, see 
Massey |2003l . An handful of field RGB stars (at F450W> 24.0 
and F450W-F814W^ 2.0) is the only population identified in 
this inner annulus which is clearly not associated with the clus- 
ter. 

The upper right and lower left panels of Fig. [9] shows the 
CMD of outer annuli of the PC field with the same area as the 
10 px < r < 160 px annulus. Even if these fields still contain 
some cluster members, their stellar mix should be fairly repre- 
sentative of the surrounding field population (compare with the 
WF2 CMD shown in the lower right panel). The comparison of 
the innermost annulus with the outer two of the same area shows 
that the supergiant population is characteristic of the cluster and 
is much less frequent in the field, suggesting an older average 
age of the field population with respect to the cluster. The com- 
parison between the morphologies of the MS is consistent this 
view. The lower right panel of the figure shows the CMD of a 
WF2 field whose area is 32 times that of the annuli described 
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Fig. 9. CMDs of different circular annuli around the center of VdBO in the PC field (see Fig. [2] above), all having the same area, 
(upper panels and lower left panel) and of the whole WF2 field, whose area is 32 times that of the PC annuli (lower right panel). The 
thin line is a Z = Zq isochrone of age 25 Myr; the heavy lines at F450W < 24.0 are 12 Gyr old isochrones of metallicity Z = 6x 10~ 4 
and Z = 6 x 10~ 3 , from blue to red, respectively. The additional isochrone plotted in the lower right panel has Z = 0.008 and age 
125 Myr. All the isochrones are from G02. 



above. The larger sampled area provides a clearer picture of the 
population mix of the M3 1 disk in the surroundings of VdBO. 
While MS and evolved stars of age (mass) similar to that encoun- 
tered in the cluster are present, the majority of the stars seem to 
have ages greater than 100 Myr. In particular the evolved stars at 
FA5QW - FUAW > 2.0 and FA5QW < 24.0 that are well fitted 
by the over-plotted 125 Myr, Z = 0.008 isochrone are not seen 
in the 10 px < r < 160 px annulus. 



The CMD of the cluster (innermost annulus) is very similar 
to that of rich Large Magellanic Cloud clusters of age ~ 30-50 
Myr, such as NGC 1711 ( Sagar et al. 11991) an d, in particular, 
NGC 1850 (Vallenari et al. [19951 Gilmozzi et al. [19951 

3.2. Supergiant Stars 

The analysis illustrated in Fig. [10] and reported in Table [3] quan- 
titatively demonstrates the presence of a significant overabun- 
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Fig. 10. CMDs of different annuli around the center of VdBO in the PC field (upper panels) and of a large area in the WF2 field 
(lower left panel), expected to sample the surrounding "field" population. An isochrone of Z = Zq and age 25 Myr is superposed 
on the upper left CMDs, as a reference. The Cr = 0.90 line is reported and a raster of labeled boxes is also over-plotted. The lower 
right panel reports the background-subtracted star counts (see Tab. [3]l in the various boxes, in units of <x, for the inner (r < 160 px, 
filled circles) and outer (160 px < r < 330 px, open circles) annuli. Zero, three and five cr levels are marked by dashed horizontal 
lines. 



dance of supergiants in the cluster with respect to the surround- 
ing field. We counted stars in the different boxes on the CMDs 
shown in Fig. [TO] sampling the upper MS (box A) and super- 
giants of blue (B), intermediate (C) and red (D) colors. The 
counts obtained in the r < 160 px and 160 px< r < 330 px 
annuli are compared with those expected from the field popula- 
tion, computed by rescaling the observed counts in the WF2 field 
by the ratio of the sampled areas. The lower right panel shows 
that in the r < 160 px annulus a clear excess of stars is present in 
all of the boxes considered. The excess of bright MS stars is very 
significant and the excess of RSGs is above the 3<x level. Even if 
the low number of stars prevents the detection of significant ex- 
cesses, the 160 px< r < 330 px annulus shows some excess with 



respect to the field in all of the considered boxes, in agreement 
with the results of Fig. [8] 

The total background-subtracted number of RSGs at- 
tributable to VdBO is =i 18. The true number is likely larger 
than this, as some RSGs are likely to reside in the innermost 
r < 10 px, which are not included in the present analysis as they 
are not well resolved in our images. According to Figer (2008) 
a richer harvest of RSGs is observed in only one known YMC 
of the Milky Way, RSGC2, with twenty-six RSG stars. RSGC1 
has fourteen, while other young clusters listed by Figer have less 
than five. RSGC2 is reported to have an age < 21 Myr, RSGC1 
has age < 14 Myr, and all the other clusters listed by Figer have 
ages < 7 Myr, i.e. younger than VdBO (see below). As noted 
above, some rich clusters of similar age are known in the LMC 
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Table 3. Star counts in the CMD boxes defined in Fig. [10] Box A samples the upper MS, boxes B, C, and D samples SG stars of blue, 
intermediate and red colors, respectively. N exp is the number of stars expected in a given box from the field population, computed 
by rescaling the observed counts in the WF2 field by the ratio of the sampled areas. The ratio between the area of the considered 
field (annulus) and the area of the WF2 field (used as representative of the field population) is reported in the last column. 
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(Vallenari et al. lT994l Brocato et al. l200TT l, but even there RSGs 
are not present in large numbers. 

3.3. Age and metallicity 

Having fixed the amount of reddening and the distance modu- 
lus to the cluster, we obtain an age estimate and an indication of 
the metallicity by comparison with isochrones from the G02 set, 
following the approach used by WH01. In Fig. [TT]we present a 
comparison with isochrones of various metallicities in the range 
|Zq ^ Z $ 2Zq. In all the panels, the isochrone that is judged 
(by eye) to provide the best-fit to the observed CMD is plotted 
as a continuous line. Dashed lines correspond to isochrones pro- 
viding strong upper and lower limits to the age estimates, which 
serve as conservative estimates of the associated uncertainties. 

The first very basic conclusion to be drawn from the reported 
upper/lower limits, is that, independent of the adopted metallic- 
ity, the age of VdBO must be within the relatively narrow range 
from 12 to 63 Myr. 

The wide range in magnitude covered by supergiant stars 
strongly indicates the presence of a wide blue loop (Massey 
2003). The super-solar isochrones clearly lack this feature, hence 
can likely be excluded as a possible solution. The larger range of 
color and magnitude covered by the Z = Zq isochrone in the 
blue loop phase seems to provide a slightly better description of 
the CMD, compared to the Z = 0.008 case. We produced a set of 
synthetic CMDs for populations having Z = 0.008, 0.019, 0.04, 
age 30 Myr and 50 Myr, Kroupa (12001 V Salpeter (1 19551 1 and 
N(m) oc m -1 ' 35 Initial Mass Function^ (IMF), using the ded- 
icated Web Tool provided by the BASTI team. After applying 
the appropriate distance modulus and reddening correction and 
transforming to the HST VEGAMAG system using the trans- 
formations by Dolphin ( 2000b|, we computed a Blue to Red 
Supergiant ratio defined as the ratio of stars having F%\4W < 
20.0 and F450W-F814W < 2.0 (B) or F450W - FS14W > 2.0 
(R). Independent of age and IMF, all the Z = 0.008 models 
have B/R < 0.26 (B/R < 0.02 mag for age = 30 Myr), while 
the observed number is B/R = 0.60 ± 0.27. The Z = 0.04 
models have 0.15 < B/R < 0.52, while the solar models have 
0.61 < B/R < 1.17. Therefore, the color distribution of SGs 
provides further quantitative support to the conclusion that the 
metallicity of VdBO is nearly solar. Adopting Z = Zq as our 
best estimate for the cluster metallicity, the age may be more 
quantitatively constrained by the comparison of the observed 
MS Luminosity Function with those predicted by models of var- 
ious ages. Fig. [12] clearly shows that an age=25 Myr model 
provides the best-fit to the observed drop in the star counts at 
F814W ^ 21.0. The result is well reproduced also if a Kroupa 
IMF is adopted. 
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Fig. 12. Comparison of the observed LF with theoretical mod- 
els from the G02 suite. Upper left panel: CMD of VDBO with 
overplotted the box adopted to select the sample of stars to be 
included in the LF. The considered radial range avoids the inner- 
most region where the completeness displays significant radial 
variations in the range of magnitudes considered. Upper right 
panel: completeness as a function of magnitude for the color and 
radial range considered. Lower panel: the observed LF (before 
completeness correction = histogram; corrected for complete- 
ness = filled circles with error bars ) is compared with models 
of different ages. Note the good fit of the drop at F8 14W - 21 .0 
achieved by the age=25 Myr model. The theoretical LF have 
been arbitrarily normalized to best match the three faintest ob- 
served points . 



Our age estimate is not expected to depend critically on the 
set of theoretical models adopted. In their thorough comparison, 
Gallart, Zoccali & Aparicio (2005) showed that there is reason- 
ably good agreement between all the theoretical isochrones they 
considered in this range of ages (i.e. < 100 Myr), if stellar mod- 
els with core overshooting are assumed. Our own (limited) set of 
experiments with Pietrinferni et al. d20041 l models also supports 
this conclusion. A few tests with a set of isochrones adopting 
the canonical treatment of convection (from Pietrinferni et al. 
2004) has shown that the adoption of such models would lead to 
younger age estimates, by a factor of ~ 4, compared to models 
including overshooting. 



9 Salpeter's IMF has N(m) oc m 235 ; Kroupa's IMF has N(m) oc m 23 
for M > 0.5M Q , and N(m) oc m' 13 for M < 0.5Mq 



Given all the above, we adopt Z — Zq as our best guess for 
the cluster metallicity, and 25 Myr as our best estimate of its 
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Fig. 11. Age estimates for VdBO for different assumptions about the total metallicity (Z). Isochrones from the G02 set are compared 
to the CMD of the cluster (10 px < r < 160 px). The best-fit isochrone is plotted as a thick continuous line while the dashed 
isochrones bracket the upper and lower limits on age. The ages and metallicities of the adopted isochrones are reported in each 
panel. The dotted lines mark the limiting magnitude as a function of color: the diagonal plume of stars just above the lines (with 
F450W-F814W> 1.5) is populated by likely RGB and AGB field stars, not associated with the cluster. 



age (see Table 5}. The mass of the stars at the TO of the best-fit 
isochrone is M TO = 9.7 Mq. 

This relatively rough age estimate is sufficient for our pur- 
poses. Our final aim is to place the cluster into a log(Age) ver- 
sus absolute integrated magnitude diagram such as that shown 
in Fig. [14] below (see also Bellazzini et al. 2008, hereafter B08, 
and references therein), to compare its stellar mass with that of 
Galactic open clusters of similar ages. The uncertainties reported 
here as the adopted upper and lower limits to the age estimates 
correspond to £ +0.3 dex in log(Age). These imply relatively 
small changes in the final estimate of the total stellar mass (a 
factor of ^ 2); the mass estimate also depends relatively weakly 
on the assumed IMF - see below - and very weakly on the metal- 
licity, at least in the range considered here, see B08). 



3.4. Integrated photometry, surface brightness profile and 
structural parameters 

Surface-brightness profile-fitting was carried out using meth- 
ods similar to those of Barmby et al. (2007). A more detailed 
description and the results of profile-fitting for the full clus- 
ter sample will be presented in Barmby et al. (2009, in prep.). 
Briefly, the two PC images in each filter were combined with the 
STScI Multidrizzle software. Intensity profiles were measured 
using the ellipse fitting routine in IRAF, on logarithmically- 
spaced isophotes centered on the intensity peaks of the clus- 
ters. The isophotal profiles were 'circularized' by converting 
the semi-major axes a of the ellipses to effective radii R e s = 
Vfl(l - e), converted to electrons s _1 arcsec~ 2 by multiplying by 
(lpixel/0.0455") 2 = 483.033 and then to intensity in L© pc~ 2 by 
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multiplying by 14.276 and 6.746 for F450W and F8 14W, respec- 
tivel)o The mean ellipticity and position angle obtained from 
the analysis of F450W and F814W images are very similar. For 
this reason we take their average as our best values, e = 0.44 
and PA=45.5°, measured from North toward East. The available 
prescription for correcting WFPC2 photometry for CTE effects 
deals only with photometry of point sources, not semi-resolved 
objects such as extragalactic star clusters; accordingly, no CTE 
corrections were made to the profiles. 

Cluster structural models were fit to the profile using the 
methods described in McLaughlin et al. ([2008 ). Before fitting 
to the data, the models were convolved with a PSF profile de- 
rived from ellipse measurements of Tiny Tim model Point Spread 
Functions (PSFs) for the center of the PC camera. We consid- 
ered the same three models used in Barmby et al. (2007): King 
( fl966l >, Wilson ( [19751 and Sersic ( fl968l . The background level 
(i.e., the intensity of the largest isophotes) was allowed to vary 
in the fitting. Fig.[T"3lshows the profile data and the best-fit mod- 
els in the two filters. Small scale bumps in the observed profile 
are likely due to individual bright stars (SGs). For the F450W 
filter the Sersic model with index n = 4.0 was the best fit. This 
model has central intensity Iq = 7.9 x 10 5 Z<j pc~ 2 and scale 
radius ro = 6.1 x 10~ 4 pc. The projected half-light radius is 
r/, = 9.12 pc (2'.'40) and total luminosity (corrected for extinc- 
tion) 1.5 X 1O 6 L . F or the F814W image, the best-fit model 
was a Wilson (1975 ) model with Wo = 11.2, central intensity 
Io = 5.0 x 1O 5 L0 pc~ 2 and scale radius ro = 0.072 pc. The pro- 
jected half-light radius is r/, = 5.60 pc (l'/47) and total luminos- 
ity 5.7 x 10 5 Lq. In the following analysis, we adopt the average 
of the two half-light radii, r h = 7.4 ± 2.5 pc (1'.'94 + 0766; the 
reported uncertainty is the standard deviation of the two values). 
It is also interesting to note that the half-light radius we have 
derived for VdBO is larger than those for the clusters listed by 
Figer (see Davies et al. 2008, r/, ^ 0.2 - 3pc), but smaller than 
NGC 1850 (r h 13 pc) and very similar to NGC 1711 (r h ^6 
pc), for exampl^iJ A summary of the adopted structural param- 
eters of VdBO is reported in Table [4] 

The derived values of the total luminosity correspond to 
M450W — -10.13 and M%iw — -10.25, respectively. Using 
Eq. 12 of Dolphin d2000bl ) these VEGAMAG magnitudes can 
be transformed to standard B and / using the appropriate co- 
efficients from his Table 7. The integrated (B - V)o color re- 
quired for the transformation has been taken from the RBC 
((B-V) = 0.05, T ab.Q] a bove), while we adopted (V-I) = 0.40 
from Maraston's (2005) model for a solar metallicity Simple 
Stellar Population (SSlF 2 !) with age of 25 Myr, as an obser- 
vational estimate of the / magnitude of VdBO was not avail- 
able (but see below). My = -9.9 is obtained from M&mw and 
My = -10.2 from M450W; we adopt the average (in flux) of the 
two, My = -10.06. This value is in excellent agreement with the 



10 This conversion assumes DN zeropoints of Z450 = 21.884, Zgi4 = 
21.528, a gain of 7 electrons DN~', and Mq F4J0W = 5.31 and 
A/0.F814W =4.14. 

11 The surface brightness profiles of these and other LMC clusters 
have been studied by Mackey & Gilmore (2003 1 who provide the pa- 
rameters of the EFF87 models that best fit the observed profiles. To de- 
rive the reported half-light radii we searched for the King 1962 model 
providing the best match to the EFF87 best-fit profile found by Mackey 
& Gilmore (2003 1, and adopted the corresponding r/,. 

12 A Simple Stellar Population is a population of stars all having the 
same age and chemical composition and having individual masses ex- 
tracted from a given Initial Mass Function (IMF); this is a practical ide- 
alized model that is generally believed to be a reasonable approximation 
of a star cluster, see Renzini & Fusi Pecci ll988l 
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Fig. 13. Intensity profiles from surface photometry in circular 
annuli from the F814W image (upper panel) and for the F450W 
image (lower panel). The continuous lines are the respective 
best-fit models, convolved with the instrumental PSF and with 
a constant background level added. For the parameters of the 
best-fit models see text. 

value of My = -10.03 listed in the RBC, and coming, in turn, 
from the photometry by Sharov et al. ( I1995I ). 

There are, however, compelling reasons to consider the es- 
timate of My obtained from our HST images as significantly 
uncertain because of the unfortunate combination of a very ex- 
tended cluster and of a very low intrinsic background level (just 
1 to 2 DN in the background sky in the original raw WFPC2 
images, particularly for the F450W filter). This guarantees that 
photometry within very large apertures will have a large uncer- 
tainty, and the resulting integrated brightness may depend on the 
details of how the code handles the background estimate in this 
photon-starved regime. 

For this reason we prefer to rely on the excellent ground- 
based material that is publicly available to obtain a reliable esti- 
mate of the total luminosity of the cluster. Existing ground-based 
photometry of VdBO taken from Sharov et al. ( 1 19951 ) is compiled 
in the RBC. However, it is possible that it was obtained adopting 
apertures that were not large enough to include the whole light 
distribution of this particularly extended cluster (see Fig. [8] and 
[L3l . We have therefore used two independent and well calibrated 
publicly available imaging surveys covering M31 to determine 
the integrated brightness of the cluster VdBO, that of Massey 
et al. (120061 hereafter M06) and the Sloan Digital Sky Survey 
(SDSS). In both cases we use an aperture with r=14'.'4. From 
the BVRI images of the former we obtained B = 14.94 ± 0.09, 
V= 14.67 ± 0.05, R = 14.45 + 0.11 and/ = 14.01 ±0.1 l B The 
SDSS - Data Release 6 (DR6, Adelman-McCarfhy et al. 120081 
g, r, and i images yielded B = 14.92, V = 14.63, R = 1 4.45, and 
/ = 14.03 using the color transformations of Lupton (2003J), in 
excellent agreement with those inferred from the M06 images. 
This is =s 0.4 - 0.6 mag brighter than those reported in the RBC. 

13 We note that these values imply (V - I)o = 0.41, adopting the red- 
dening law by Dean, Warren & Cousins dl978t . in excellent agreement 
with the prediction, used above, of (V - /)o from Maraston's (2005) 
model for a solar metallicity SSP of age 25 Myr. 
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Fig. 14. Integrated V mag and total mass as a function of age for 
various clusters. Galactic Open Clusters (OC, from the WEBDA 
database) are plotted filled circles, Galactic Globular Clusters 
(GC, My from the most recent version of the Harris (1996) cat- 
alogue, i.e. that of February 2003; the ages have been arbitrarily 
assumed to be 12.0 Gyr for all the clusters) are plotted as x sym- 
bols. VdBO is represented as a crossed square at My = -10.42, 
from Tab. [4] The continuous lines are fixed-stellar-mass models 
from the set by Maraston (1998, 2005) for SSPs of solar metal- 
licity, with a Salpeter's Initial Mass Function (IMF) and inter- 
mediate Horizontal Branch morphology. Note that in this plane, 
the dependence of the models from the assumed IMF, metallicity 
and HB morphology is quite small (see B08). The outlier OC at 
log Age^ 9.0 is Tombaugh 1. The long dashed line is the VDBO 
evolutionary track including the mass loss by dynamical effects 
according to the formulas by LG06. The cluster is expected to 
dissolve within < 4 Gyr from the present epoch. 

In Sect. 4 we will show that the J,H,K magnitudes of VdBO also 
become brighter by ~ 0.2 - 0.5 mag after increasing the adopted 
aperture from r =5'.'0 to 15'.'0. 

Given all the above, we adopt the r -14'.' 4 aperture photom- 
etry measured on M06 images as our preferred values, reported 
in Table|H below. In particular V = 14.67 + 0.05 is our final best 
estimate of the integrated V magnitude of VdBO, correspond- 
ing to My = -10.42 ± 0.20; these values will be adopted in the 
following analysis. 

4. Summary and discussion 

We have outlined the data reduction and scientific analysis strat- 
egy that we adopt for our HST-WFPC2 survey of M3 1 candidate 
YMCs, whose complete results will be presented in future con- 
tributions. As an exemplary case, we have described the study of 
the cluster VdBO. We have found that VdBO is a very bright and 
extended cluster of approximately solar metallicity and of age 
~ 25 Myr, with a rich population of blue and red supergiants. 

Having clearly ascertained that VdBO is a real cluster, it re- 
mains to be established if it is more similar to ordinary open clus- 



Fig. 15. The same as Fig. [14] but for near infrared colors. 
Integrated magnitudes of GCs are taken from Cohen et al. 
( 120071) ; the IR magnitudes for VdBO are taken from Tab. 4. The 
dotted lines are M = 10 4 M© and M = 10 5 M o iso-mass models 
assuming a Kroupa 2001 IMF instead of a Salpeter ( 119551 ) IMF, 
plotted here to illustrate the weak effect of assumptions on IMFs. 



ters of the Milky Way than to to the Young Massive Clusters that 
may be considered as possible precursors of "disk globulars". 
The similarity with LMC objects typically classified as "Young 
Globular Clusters" such as NGC1850 (see Sect. 3., above) is 
quite remarkable and it suggests that VdBO is not an ordinary 
OC (but see also point 1, below). 

A more general way to compare clusters of different ages, 
taking into account the fading of the luminosity of SSPs as they 
age, it is to plot them into a diagram comparing age to some 
indicator of the stellar mass of the cluster (see, for example, 
Whitmore, Chandar & Fall 120071 Gieles, Lamers & Portegies- 
Zwart 2007 and de Grijs, Goodwin & Kouwenhoven 20081 for 
recent applications and references). Here we adopt log(Age) vs. 
absolute integrated magnitude as in B08. 

In Fig. [14] VdBO is compared with Galactic Open Clusters 
(data taken from the WEBDA databasfl), with Galactic 
Globular Clusters (from the latest version of Harris 1996 assum- 
ing a uniform age of 12 Gyr, a reasonable approximation for our 
purpose), and with a grid of SSP models with solar metallicity 
and Salpeter's IMF from the set by MarastorH ( Tl998ll2003T l. As 
a SSP ages massive stars die while the mass of the most luminous 
stars decreases (passive evolution). Keeping the total mass fixed, 
the luminosity of the population fades and, as a consequence, the 
stellar mass-to-light (M/L) ratio increases. The continuous lines 
plotted in Fig. [14] describe the passive evolution of SSPs of var- 
ious (stellar) masses: under the adopted assumptions the mass 
of a cluster of given age and My can be read from the grid of 
iso-mass tracks. 



http : //www . univie . ac . at/webda/integre . html 
|http : //www-astro . physi cs . ox . ac . uk~ maraston/ 
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The path of the track passing through the cluster shows what 
its luminosity will be in the future if the cluster did not lose stars 
through dynamical processes (evaporation, tides, ecc.). The lat- 
ter is clearly not the case in general, and in particular for VdBO. 
In addition to the relatively mild evaporation driven by two body 
encounters, it will suffer from the strain of the M3 1 tidal field 
and from encounters with Giant Molecular Clouds (GMC), as 
the cluster is embedded in the dense thin disk of M3 1 (Lamers 
& Gieles 2006 hereafter LG06, and references therein). To take 
these effects into account we used the analytical approach pre- 
sented by LG06 to produce an evolutionary track including the 
cluster mass loss by stellar evolution, galactic tidal field, spi- 
ral arm shocking, and encounters with giant molecular clouds, 
plotted in Fig. [14] as a long-dashed curve. The LG06 formulas 
describe the evolution of a cluster located within the Milky Way 
(thin) disk at the Solar circle. They should provide a reasonable 
approximation for VdBO which lies in the disk of M31, at a sim- 
ilar distance from the center of a similarly massive spiral galaxy 
(van den Bergh 2000). The required inputs are the cluster mass, 
for which we adopted the value that can be read from the SSP 
grid of Fig. [14] (see below), and the half-light radius, which we 
obtained in Sect. 3.3, above (see Tab.[4]l. The initial expulsion of 
gas not used in star formation may lead young clusters (age < 50 
Myr) to lose their virial equilibrium and it may represent an ad- 
ditional relevant factor driving toward the destruction of clusters 
like VdBO that is not included in the LG06 approach (Bastian & 
Goodwin 120061 Goodwin & Bastian 120061 Bastian et al. [2008l . 

Fig.Q3]is worth of some detailed considerations: 

1. Independently of the exact value of My adopted, VdBO is 
significantly brighter (<: 1 mag) than Galactic OCs of simi- 
lar ages, actually it is brighter than Galactic OCs of any age. 
The same is true also if all other known M3 1 OCs are consid- 
ered (Hodge [1979l Krienke & Hodge 120071 120081 . However 
it should be noted that the population of disk clusters in M3 1 
may be so huge (~ 80000 clusters, according to Krienke & 
Hodge 2007) that even the extreme tails of the luminosity 
distribution may be populated. (This should not be the case 
for the LMC, for example, as it is orders of magnitude less 
massive than M31). Hence it is premature to draw a conclu- 
sion from an individual cluster; when the whole sample is 
analyzed we will get a deeper insight on the actual nature of 
VdBO. 

2. Assuming the RBC value for the integrated V magnitude, 
E(B-V)=0.0 instead of E(B-V)=0.2 and a grid of iso-mass 
tracks adopting a Kroupa IMF, we can obtain an extremely 
conservative strong lower limit to the stellar mass of VdBO, 
M = 2.4x 1O 4 M0. Under the same assumptions but adopting 
the best-fit value E(B-V)=0.2 we obtain M = 6.5 x 10 4 M© 
with a Salpeter IMF and M = 4.2 x 1O 4 M0 with a Kroupa 
IMF. These are at the threshold between the OC and GC 
mass distributions (see van den Bergh & Lafontaine 1 19841 
and B08) and also at the upper end of the mass distribution of 
Galactic YMC (see Figer 2008 and Fig. [16] below). The con- 
clusion that VdBO is much more massive than MW clusters 
of similar ages seems inescapable, unless extreme IMFs are 
considered (i.e. IMF truncated at low masses, see Sternberg 
[19981 . 

3. If My = -10.42 is adopted, as obtained from large aper- 
ture ground-based V photometry in Sect. 3.4, the total stel- 
lar mass is M = 9.5 x 10 4 M Q with a Salpeter IMF and 
M = 6.0 x 10 4 M o with a Kroupa IMF. 

4. The evolutionary tracks including the LG06 treatment of 
mass-loss by dynamical effects show that, independent of 



the actual mass (within the range outlined above), it is un- 
likely that the cluster VdBO would survive for an Hubble 
time. Hence it is very probable that it will never have the op- 
portunity to evolve into a classical (faint) GC. The disruption 
timescale is dominated by encounters with GMCs; consider- 
ing this effect alone (Eq. 7 of LG06) the cluster is predicted 
to dissolve within ^3.6 Gyr if its mass is M = 9.5 x 10 4 M o , 
as obtained from our best estimate of the integrated V mag- 
nitude and assuming a Salpeter's IMF. 
5. In the same grid of Fig. [14]and under the same assumptions 
the masses of the BLCCs observed by WH01 - adopting their 
age estimates - range from 8.0 x 1O 3 M0, (G293) in the realm 
of OCs, to 2 x 10 4 M o (G44 and G94) and 8 x 10 4 M o 
(G38), very similar to that of VdBO and significantly larger 
than OCs of similar ages. 

To obtain independent and more robust estimates of the 
present-day stellar mass of VdBO we used the Near Infrared 
(NIR) version of the log Age vs. absolute integrated magnitude 
plane. In Fig. [15] J,H and K absolute magnitudes of VdBO ex- 
tracted from the Extended Sources Catalogue (XSC) of 2MASS 
are compared with Maraston's SSP models of solar metallicity 
and Salpeter's (continuous lines) or Kroupa's (dotted lines) IMFs 
and with Galactic GCs (from Cohen et al. 2007 , ages assumed 
as above NIR integrated magnitudes for significant samples 
of OCs are not available, at present. To account for the whole 
extent of the cluster we extracted r = 15" aperture photometry, 
that is provided in the XSC, instead of the r = 5" adopted in the 
RBC, see Tab. Hand Tab.©. 

NIR magnitudes are more reliable mass tracers than visual 
magnitudes as NIR M/L ratios are smaller and have smaller vari- 
ations with age, compared to optical M/L ratios. For example, 
according to Maraston (1998] 120051) models, a solar metallic- 
ity Salpeter-IMF SSP at Age = 10 Gyr has (M/L) v =5.5, while 
(M/Lk=1.4; the same SSP has * 0.55 while * 

0.13. The independent estimates of the stellar mass from J,H, 
and K magnitudes are essentially identical, ranging from 6 to 9 
xlO 4 M0, assuming a Salpeter IMF, and from 4 to 5.5 x10 4 Mq, 
assuming a Kroupa IMF. These estimates are in fair agreement 
with those obtained from the integrated V photometry. 

Finally, in Fig. [16] we compare VdBO with Galactic OCs, 
GCs and YMC, plus the BLCCs studied by WH01, in the log of 
the stellar mass versus log of the half-light radius plane (similar 
to Mackey & van den Bergh [2005] and Federici et al. 120071 . The 
radii at which the break in the profile (core/corona transition) of 
Galactic OCs (from Kharchenko et al. 2005} occurs is taken as 
a strong upper limit for their r/,. VdBO has a typical size that is 
larger than both OCs and YMCs, and is similar to that of several 
MW GCs of comparable mass. 

In conclusion, we can say that VdBO seems a remarkable 
cluster in several of its properties when compared to the other 
known disk clusters of the Milky Way and M31. In this paper 
we have presented the data reduction, data analysis and diagnos- 
tics that will be applied to the whole survey sample and that will 
allow us to put VdBO and the other clusters in the more gen- 
eral context of the star cluster populations in the disk of spiral 
galaxies. 
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through the PRIN 2007 grant CRA 1.06.10.04 "The local route to galaxy forma- 
tion...". P.B. acknowledges research support through a Discovery Grant from the 



16 For J,H,K colors we adopt Aj = 0.871£(5 - V), A H = .540£(fi - 
V), and A K = 0.346£(B - V), from Rieke & Lebof sky 19851 The J,H,K 
absolute magnitudes of the Sun are taken from Holmberg, Flynn & 
Portinari l2006l 
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Table 4. Newly derived coordinates, half-light radius, integrated 
magnitudes, reddening, age and metallicity for the cluster VdBO. 
The origin of each parameter is described in the last column. 



par 


value 


note 


<?/2000 


00* 4(y» 29.4 s 


from 2MASS-XSC 


<5/2000 


+40° 36' 15.2" 


from 2MASS-XSC 


ri, 


1'.'93 ± a.' 66 


from intensity profile (i.p.) fit 


B 


14.94± 0.09 


r=14'.'4 ap. phot, on M06 images 


V 


14.67± 0.05 


r=14'.'4 ap. phot, on M06 images 


R 


14.45+ 0.11 


r=14'.'4 ap. phot, on M06 images 


I 


14.01+ 0.11 


r=14'.'4 ap. phot, on M06 images 


J 


13.26+ 0.07 


r=15'.'0 ap. phot, from 2MASS-XSC 


H 


12.76+ 0.12 


r=15'.'0 ap. phot, from 2MASS-XSC 


K 


12.77+ 0.15 


r=15'.'0 ap. phot, from 2MASS-XSC 


age 


25 Myr 


value of adopted best-fit isochrone 


Z 


0.019 


value of adopted best-fit isochrone 


E(B-V) 


0.20 


adopted best-fit value 



Fig. 16. VdBO (crossed square) is compared to other clusters 
in the logarithm of the mass vs. logarithm of the half-light- 
radius plane. Filled circles are Galactic GCs from Mackey & 
van den Bergh (2005). Arrows are Galactic OCs: we plot the 
radii where a break in the surface brightness profile occurs, taken 
from Kharchenko et al. (2005, their "core radii"). These should 
be considered as upper limits for actual r/,, which are not avail- 
able for most OCs. The masses of the OCs have been com- 
puted using the grid of SSP models shown in Fig.[l4l while for 
GCs we adopted age=12.0 Gyr and a grid of SSP models hav- 
ing [Z/H] = -1.35. Open pentagons are the clusters studied by 
WH01 . Open triangles are the massive young MW clusters listed 
by Figer (2008); masses and radii are taken from his Table 1. 
Note that the radii reported by Figer for these clusters are not 
half-light-radii, however they should be a reasonable proxy. The 
good match between the two quantities has been verified in the 
case of RSG1, for which Figer report r = 1.3 pc, and Davies et 
al. d2008b obtain r h = 1.5 ± 0.3 pc. 
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